A feature-based mesh adaptation for the unsteady high speed compressible flows in complex three-dimensional domains

We propose an unstructured mesh adaptation approach for unsteady high speed compressible Navier–Stokes applications involving blasts and explosions with the presence of strong shock waves propagating in three dimensional complex domains. The idea is to identify the locations of critical physics locally and then re-mesh these regions based on solution derived metrics. The approach ensures both geometry fidelity and mesh validity, especially for areas near complex geometries, a task that is always a challenge in mesh adaptation. The proposed adaptivity is applied for simulations of blast wave propagations and compared with available data in literature. The results show that the proposed method is fully robust and efficient for computational fluid dynamics (CFD) problems in complex three-dimensional domains

EM modelling of arbitrary shaped dispersive chiral dielectric objects using a 3D leapfrog scheme on unstructured meshes

The standard Yee FDTD algorithm is widely used in computational electromagnetics because of its simplicity and divergence free nature. A generalization of this classical scheme to 3D unstructured co-volume meshes is adopted, based on the use of a Delaunay primal mesh and its high quality Voronoi dual. This circumvents the problem of accuracy losses, which are normally associated with the use of a staircased representation of curved material interfaces in the standard Yee scheme. The procedure has been successfully employed for modelling problems involving both isotropic and anisotropic lossy materials. Here, we consider the novel extension of this approach to allow for the challenging modelling of chiral materials, where the material parameters are frequency dependent. To adequately model the dispersive behaviour, the Z-transform is employed, using second order Padé approximations to maintain the accuracy of the basic scheme. To validate the implementation, the numerical results produced are compared with available analytical solutions. The stability of the chiral algorithm is also studied

A 3D Unstructured Mesh FDTD Scheme for EM Modelling

The Yee finite difference time domain (FDTD) algorithm is widely used in computational electromagnetics because of its simplicity, low computational costs and divergence free nature. The standard method uses a pair of staggered orthogonal cartesian meshes. However, accuracy losses result when it is used for modelling electromagnetic interactions with objects of arbitrary shape, because of the staircased representation of curved interfaces. For the solution of such problems, we generalise the approach and adopt an unstructured mesh FDTD method. This co-volume method is based upon the use of a Delaunay primal mesh and its high quality Voronoi dual. Computational efficiency is improved by employing a hybrid primal mesh, consisting of tetrahedral elements in the vicinity of curved interfaces and hexahedral elements elsewhere. Difficulties associated with ensuring the necessary quality of the generated meshes will be discussed. The power of the proposed solution approach is demonstrated by considering a range of scattering and/or transmission problems involving perfect electric conductors and isotropic lossy, anisotropic lossy and isotropic frequency dependent chiral materials

Development of an axisymmetric parallel solution algorithm for membrane separation process

A novel parallel technique that couples the lattice–Boltzmann method and a finite volume scheme for the prediction of concentration polarisation and pore blocking in axisymmetric cross–flow membrane separation process is presented. The model uses the Lattice–Boltzmann method to solve the incompressible Navier–Stokes equations for hydrodynamics and the finite volume method to solve the convection–diffusion equation for solute particles.Concentration polarisation is modelled for micro–particles by having the diffusion coefficient defined as a function of particle concentration and shear rate. The model considers the effect of an incompressible cake formation. Pore blocking phenomenon is predicted for filtration membrane fouling by using the rate of particles arriving at the membrane surface.The simulation code is parallelised in two ways. Compute Unified Device Archi- tecture(CUDA) is used for a cluster of graphical processing units(GPUs) and Message Passing Interface(MPI) is utilised for a cluster of central processing units(CPUs), with various parallelisation techniques to optimise memory usage for higher performance. The proposed model is validated by comparing to analytical solutions and experimental result

A novel implementation of computational aerodynamic shape optimisation using Modified Cuckoo Search

This paper outlines a new computational aerodynamic design optimisation algorithm using a novel method of parameterising a computational mesh using `control nodes'. The shape boundary movement as well as the mesh movement is coupled to the movement of user--defined control nodes via a Delaunay Graph Mapping technique. A Modified Cuckoo Search algorithm is employed for optimisation within the prescribed design space defined by the allowed range of control node displacement. A finite volume compressible Navier--Stokes solver is used for aerodynamic modelling to predict aerodynamic design `fitness'. The resulting coupled algorithm is applied to a range of test cases in two dimensions including aerofoil lift--drag ratio optimisation intake duct optimisation under subsonic, transonic and supersonic flow conditions. The discrete (mesh--based) optimisation approach presented is demonstrated to be effective in terms of its generalised applicability and intuitiveness

Advances in co-volume mesh generation and mesh optimisation techniques

This paper introduces developments in modified techniques for the generation of unstructured, non-uniform, dual orthogonal meshes which are suitable for use with co-volume solution schemes. Two new mesh generation techniques, a modified advancing front technique and an octree-Delaunay algorithm, are coupled with a mesh optimisation algorithm. When using a Delaunay–Voronoi dual, to construct mutually orthogonal meshes for co-volume schemes, it is essential to minimise the number of Delaunay elements which do not contain their Voronoi vertex. These new techniques provide an improvement over previous approaches, as they produce meshes in which the number of elements that do not contain their Voronoi vertex is reduced. In particular, it is found that the optimisation algorithm, which could be applied to any mesh cosmetics problem, is very effective, regardless of the quality of the initial mesh. This is illustrated by applying the proposed approach to a number of complex industrial aerospace geometries